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Abstract 

We have studied the wave properties of the cold and isothermal 
plasma in the vicinity of the Schwarzschild black hole event horizon. 
The Fourier analyzed perturbed 3+1 GRMHD equations are taken on 
the basis of which the complex dispersion relations are obtained for 
non-rotating, rotating non-magnetized and rotating magnetized back- 
grounds. The propagation and attenuation vectors along with the re- 
fractive index are obtained (shown in graphs) to study the dispersive 
properties of the medium near the event horizon. The results show 
that no information can be obtained from the Schwarzschild magneto- 
sphere. Further, the pressure ceases the existence of normal dispersion 
of waves. 

Keywords: 3+1 formalism, GRMHD equations, dispersion relations. 
PACS numbers: 95.30.Sf, 95.30. Qd, 04.30.Nk 

1 Introduction 

Black holes are solutions of the Einstein field equations (EFEs) and Scliwarzschild 
solution is the simplest black hole solution. The most physical black hole so- 
lution is the Kerr metric which is axisymmetric and rotating. This reduces 
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to the Schwarzschild black hole when its angular momentum decreases to 
zero by energy extraction [T]. Around the Schwarzschild event horizon, the 
powerful gravitational force pulls the magnetized plasma of the surrounding 
space towards the event horizon in the form of an accretion cloud. The ac- 
creting plasma creates a magnetic field. The dynamical effects of accreting 
magnetospheric plasmas near a black hole's event horizon enforce us to use 
the theory of general relativistic magnetohydrodynamics (GRMHD) to study 
them. 

The gravitomagetic waves and gravitational perturbations in the black 
hole regime have always been of interest. Many people [2]-[5] studied pertur- 
bations in the Schwarzschild regime. To make the results of GR accessible 
for asytrophysicists, Arnowitt, Deser and Misner [B] developed a formula- 
tion (ADM 3+1 formalism) which splits the four-dimensional spacetime into 
three-dimensional hypersurfaces labeled by time. This makes the results of 
GR comparable with those of Newtonian physics. Several authors [7]-[T2] 
applied this formalism to judge different aspects in GR. 

Thorne and Macdonald [13]- [SI extended the formulation to electromag- 
netic fields of black hole theory. Holcomb and Tajima [15J, Holcomb [TB] 
and Dettmann et al. [17] investigated the wave propagation in Friedmann 
universe. Buzzi et al. [IE] discussed one-dimensional radial propagation of 
transverse and longitudinal waves, in two component plasma, close to the 
Schwarzschild event horizon. The stationary symmetric GRMHD theory of 
black holes was developed by Zhang [TH]. He also investigated the behavior 
of perturbations of cold plasma in the ergosphere of Kerr black hole [20] . 
Recently, Sharif and Umber [21]- [25] have found some interesting wave prop- 
erties of cold and isothermal plasmas (with constant rest-mass density) in 
the vicinity of Schwarzschild black hole event horizon. They have evaluated 
real and complex wave numbers using 3+1 GRMHD equations. 

In this paper, we shall consider variable rest-mass density to study the 
wave properties of cold and isothermal plasmas in Schwarzschild planar ana- 
logue by using complex wave vector components. The paper has been orga- 
nized as follows: Section 1 contains the description of planar analogue with 
respective background, perturbation and Fourier analysis assumptions. Sec- 
tions 2, 3 and 4 contain the dispersion relations obtained for cold plasma liv- 
ing in non-rotating, rotating non-magnetized and rotating magnetized back- 
grounds whereas sections 5, 6 and 7 constitute the dispersion relations for 
isothermal plasma living in the same backgrounds. All this has been done 
for variable mass density and pressure. The results will be discussed in the 
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last section. 



2 Schwarzschild Planar Analogue and Rela- 
tive Assumptions 



The Schwarzschild planar analogue can be described as [20] 

ds^ = -a^{z)dt^ + dx^ + dy'^ + dz^. (2.1) 

The directions z, x and y are analogous to Schwarzschild's r, (f) and 9 re- 
spectively. In this planar analogue, we assume the existence of cold and 
isothermal plasmas with the respective equations of state 



fi=—, fi= , 2.2 

Po Po 

where po, p and p denote the rest, moving mass densities and pressure re- 
spectively. 

In non-rotating background, the perturbed flow of fluid is only along z- 
axis. The FIDO measured magnetic field and velocity are given by 

B = Be,, V = uiz)e,. (2.3) 

In rotating background, FIDO measured velocity and magnetic field are two- 
dimensional. In planar analogue, these quantities can be expressed in xz- 
plane by the following expressions 

B = B[X{z)e, + e,], V = V{z)e, + u{z)e,. (2.4) 

Here A, u and V are related to each other by the following equation 



V = + An, (2.5) 

a 

where is an integration constant. 

We assume linear (first order) perturbations in flow variables (mass den- 
sity p, pressure p, velocity V and magnetic field B) of the fluid. 

p = p^ + 6p = + pp, p = p° + 6p = p^ + pp, 

V = V° + 5V = V° + V, B = B° + 5B = B° + 5b, (2.6) 
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where p°, p, V and B are unperturbed quantities, Sp, Sp, S\ and (5B rep- 
resent perturbed quantities. The dimensionless perturbed quantities p, p, v 
and b can be written as follows 

p = p{t,z), p = p{t,z), 

h = — = b^{t,z)e^ + b^{t,z)e^. (2.7) 

For Fourier analysis, we assume the harmonic space and time dependence of 
perturbations 

p{t, z) = cie-'(^*-^^), p{t, z) = C2e-'-^^'-^'\ 
V, {t,z)= C3e-'("*-'=^) , v^{t,z) = ae-'^^'-^'^ , 

z) = C5e-'(-*-'=^), K{t, z) = C6e-^('^*-^^). (2.8) 



3 3+1 Perfect GRMHD Equations in Planar 
Analogue 

For the Schwarzschild planar analogue, the perfect GRMHD equations can 
be written as ^21j.^26j 

rfVK 

_ = V X (aV X B), (3.1) 

V.B = 0, (3.2) 

do dV 

^ + (aV).Vp + P7'V.— + p72v.(aV.V)V 

+pV.(aV) = 0, (3.3) 

{(Pf + 1^)5,- + pi'v.y, - + v.v)y^- 

= -p^^a, - p,, +-^(V X B),V.(V X B) - -i^(aB)2„ 

Arc OTTO;"^ 

+^(aB,)„ B:* - J-|B X {V X (V X (aV x B))}]„ (3.4) 
+ v.V)(.p„) - i| + 2p„,,V.(i| + V.V)V 
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+2po/x7'(V.a) + po/x7'(V.V) + ^[(V x B).(V x (aB)) 

+ (VxB).^(VxB)]. (3.5) 

The cold plasma model can be described by Eqs. fl3.1l) - fl3.4p whereas isother- 
mal plasma needs Eg. (13.51) for a better interpretation. 

In a stationary symmetric background, the observers move along sym- 
metric directions and thus do not see any change in flow around themselves. 
For such magnetized fluids, Phinney proposed the stream functions h, I and 
e il9j, i.e., 

— Attu 

h= — - — , l = 'yV — aXB/h, 6 = 70;, (3.6) 
13 

where the rest-mass conservation law in three dimensions is apo'-fu = —1 
with specific enthalpy /i = 1 and the constant = 0. 



4 Non- Rotating Background with Cold Plasma 

The Fourier analyzed perturbed GRMHD equations in non-rotating back- 
ground (Eqs.(3.15)-(3.18) of [21]) are 



C5 = 0, 

a 

ikc^ = 0, 



Cl 



a 



+ Lku] +Oi-{{l+ Yu^)ik - (1 + 7^M^)(1 - 2Yu^) 



u' 
u 



LUJ 2 

— 7 u 
a 



Ci7 [a^ + uu\l + 7 M )} + C3 



7^(1 + 7V) 



-LUO 



a 



Luk 



+-i\u'{l + 7 V) (1 + 47 V) + 2-i\a,]\ = 0. 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



Equations (14. ip and ( 14. 2p show that C5 = which means that there are no 
perturbations in magnetic field. Thus the non-magnetized as well as mag- 
netized backgrounds admit the same perturbed Fourier analyzed GRMHD 
equations. 
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4.1 Numerical Solutions 



In order to find the numerical solutions, we assume the time lapse a = 
^tanh(lOz). For stationary flow, we take 07 = 1 which implies that 7 = 
1/a. For the inflow of fluid into the black hole event horizon we take u = 
— Vl — tt^- Using these assumptions, the mass conservation law in three 
dimensions gives P = These quantities satisfy the GRMHD equations 
for the region 0.5 < 2; < 10. For these assumptions, e = 1. 

The dispersion relation from Eqs. (14.31) and (14. 4p can be obtained (using 
Mathematica) by equating determinant of the coefficients of constants ci and 
C3 to zero P?]. This determinant leads to a complex dispersion relation of 
the form 

Ai{z, uj)k'^ + ^2(2:, uj)k + y43(z, uj) + L{Ai{z, uj)k + A^{z, u;)} = (4.5) 

which gives two complex values of k. Since k is the z-component of the wave 
vector, it gives quantities corresponding to 2;-direction. The real and imag- 
inary parts of k yield the propagation and attenuation vectors respectively. 
The propagation vector gives refractive index on the basis of which the mode 
of dispersion can be found. The sinusoidal expression then takes the form 
^~,{u>t-k^z-Lk2z) ^ ^-,[^t-k^z)-k2z ^ ^Yieie ki = Re{k) and ^2 = Im(A;). The two 

values of k are shown in Figure 1 and 2. 

Figure 1 indicates that the propagation vector decreases with the increase 
in z and u. It means that the waves propagate slowly as they move away from 
the event horizon and as their angular frequency increase. The attenuation 
vector admits positive and negative values randomly which shows random 
damping and growth of waves. The refractive index is greater than one 
and its variation with respect to angular frequency is positive in the region 
1.9<z<10, l<ti;<10 which gives that the region admits normal 
dispersion of waves [28j. Figure 2 shows that the attenuation vector takes 
positive values at random points. The propagation vector decreases with the 
increase in angular frequency. The attenuation vector takes random values 
which means that the waves damp and grow randomly with the increase in 
angular frequency. The refractive index is greater than one but £ < m 
most of the region giving anomalous dispersion of waves [28] . 

In Figures 1 and 2, the propagation vector takes negative values which 
shows that the waves move towards the black hole event horizon. In addition, 
a small region near the event horizon admits decrease in refractive index with 
the increase in z. 
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Surface Plot of Propagation Vector 




Figure 1: The waves are directed towards the event horizon. Dispersion is 
normal in most of the region 



Surface Plot of Propagation Vector 




Figure 2: The waves move towards the event horizon. Most of the region 
admits anomalous dispersion of waves 
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5 Rotating, Non-magnetized Background with 
Cold Plasma 

For the rotating non-magnetized background, the Fourier analyzed perturbed 
perfect GRMHD equations are given by Eqs.(4.13)-(4.15) of 





^ —LUJ 








y a 


- Luk^ 


+ C3 



-LUJ 



a 



Yu + (1 + Yu')ik - (1 - 2Yu') 



u 



x(l +7V)- + 27V\/\/' 
u 

x{(l + 2YV^)V' + 2YuVu'}\ = 0, 
ci7^u{(l + YV^)V' + YuVu'} + C37^ 







+ C47' 


( h Lku ) 




A a / 



(5.1) 



— iCiJ 



a 



+ m/c 7 



+{(1 + 27V) (1 + 27V^) - 7V2}y' + 27^(1 + 2-i'^u^)uVu'] 



+C4 



— LUJ 



a 



+ Lku 7^(1 + 7V2) + Yu{{l + 47V2) W 



+4l^l^'(l + 7V2)}] =0, 

Ci7^{az + (1 + 7^m^)mm' + 7Vl/r'} + C3 



(5.2) 



7^(1 + 7V) 



a 



+Luk) + 7^[m'(1 + 7^M^)(1 + 47^2) + 2M72{a, + (1 + 2-f\^)VV'}]] 



Luk ] uV + u^V'{l + A-f''V 



+2V(a^ + (1 + 2-f\^)uu')] = 0. 
5.1 Numerical Solutions 



(5.3) 



We assume the same time lapse which we have assumed for non-rotating 
background. The stationary flow assumption aj = 1 with V = u leads to 

Thus the mass conservation law in three dimensions gives 



u 



p = —-. These values satisfy the perfect GRMHD equations for the range 
0.75 < z < 10. In rotating non-magnetized background, the constant values 
of u = V = —0.703562 make the flow constants I = —7.03562 and e = 1. 
Consequently, we obtain the following form of complex dispersion relation 
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Surface Plot of Propagation Vector 
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Figure 3: The waves move towards the event horizon. Most of of the region 
shows normal dispersion of waves 

cubic in k. 



This dispersion relation gives three complex values of k with corresponding 
graphs shown in Figures 3-5. 

In Figure 3, the propagation vector decreases with the increase in z. The 
attenuation vector decreases with the increase in 2; in a small region near 
the event horizon. Thus, in the small region near the event horizon, the 
waves damp as they move towards the event horizon. The refractive index 
is greater than one and ^ > except for the region 2.8<2;<10, 1< 
uj < 2.1 which shows normal dispersion of waves. Figure 4 shows that the 
attenuation vector takes random negative values with increase in angular 
frequency. In a small region near the event horizon, the attenuation vector 
increases with increase in z. This indicates that the waves grow as they move 
towards the event horizon. The refractive index is greater than one and its 
variation with respect to angular frequency is negative which shows that the 
whole region admits anomalous dispersion of waves. Figure 5 indicates that 



Ai{z)k'^ + A2{z, uj)k + A^{z, u) + L{Ai{z)k^ + A^{z, uo)k'^ 




(5.4) 
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Surface Plot of Propagation Vector 



Surface Plot of Attenuation Vector 




Figure 4: The waves move towards the event horizon. The region shows 
anomalous dispersion 




Figure 5: The waves are directed towards the the event horizon. A small 
region near the event horizon admits normal dispersion of waves 
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the attenuation vector takes negative values. Near the event horizon, the 
attenuation vector decreases with the increase in z. It means that waves 
grow near the event horizon. The refractive index is greater than one and its 
variation with respect to angular frequency is greater than zero in the region 
0.75 < z < 1.1, 1.5 X 10"^ < < 10 which admits normal dispersion of 
waves. In rest of the region, dispersion is anomalous and normal randomly 
due to variation in the refractive index with respect to angular frequency 
which attain positive and negative values randomly respectively. 

In Figures 3-5, the propagation is negative in the whole region. This 
quantity decreases with the increase in angular frequency. In Figures 4 and 
5, a small region near the event horizon shows that the refractive index 
decreases with the increase in z. 

6 Rotating Magnetized Background with Cold 
Plasma 

In this background, the Fourier analyzed form of the perturbed GRMHD 
equations (Eqs.(5.15)-(5.20) of [2lj) is given as follows 

C4(a' + ika) — csKaA)' + ikaX) + c^iaV)' — c^{{au)' — loj 

+Lkua} = 0, (6.1) 

cr,{^^ + tku) = 0, (6.2) 
a 

c^tk = 0, (6.3) 

ci(— ^ + Luk) + C3[— + (1 + l'^u^)ik - (1 - 27V) 
a a 

x(l + ^W)- + 2^^u^VV'] + c^i'^f^V + ikuV 
u a 

+u{{l + 27V2)V' + = 0, (6.4) 

Ci7V{(l + 7^2)^' + -f^uVu'} - — C6{(1 - u^)ik + -(1 - u^) 

Atx a 

-uu I + C3 -(p7 uV — ) h (p7 uV + — — jiku 

+P7^{(1 + 2-f\^){l + 27V2) - j^V^}V' + 2p7^(l + 2-f\^)uVu' 

+p^{aXy] + C4[-{pf (1 + 7V2) + + {pf (1 + I'V) 

ilia 4tc a 
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td2 

- — }iku + P7^m{(1 + 47V2)W + WV'{1 + 
47r 

B'^ua' 

= 0, 6.5 

cifn^ia, + (1 + 7V)W + 7='mVF'} + C3[-{p7'(l + 7V) 
+ ^}- + {fri\l + 7V) - ^}iku + {p7^K(l + 7V) 

x(l + 47^2) + 2u^\a, + (1 + 27\2)\/\/')} - ^(«A)'}] 

+C4[-(p7^mI^ - + (P7^?^'^ + ^^)i^ku 

An a An 

+{p7'{n V'(l + 47V2) + 2V{a, + (1 + 272^2)^)} + ^^!^}] 

47rQ; 

+ ^C6[A(1 - n')iA; + A-(l - u^) - \uu' + = 0. (6.6) 

47r a a 

Equations (16. 2p and (16. 3p give C5 = which shows that there are no pertur- 
bations in 2;-component of magnetic field. 



6.1 Numerical Solutions 

We assume the same values of time lapse, density, x and ^-component of 
velocity which we have mentioned in the previous section. Substituting u = 
V in Eq. (l2.5p with the assumption that = 0, it gives that A = 1. Also we 

take B = \J'^- These values satisfy the GRMHD equations for the region 

1 < z < 10. For this region, the numerical values of the fiow constants (due 
to constant u = V = —0.703562, B and A) are e = 1, I = —6.83562 and 
h = -2.50663. 

We obtain a complex dispersion relation quartic in k 

Ai{z)k'^ + A2{z, uj)k^ + A^{z, uj)k'^ + A^{z, uj)k + A^{z, uj) 
+L{AQ{z)k^ + Aj{z, uj)k'^ + Aii{z, uj)k + A^{z, tu)} = 0. (6.7) 

The corresponding graphs are shown in Figures 6-9. 

Figure 6 shows that the propagation vector decreases with the increase 
in UJ and z. The attenuation vector attains random values in the region 
2.1 < 2 < 10, < C(j < 10 showing that waves damp and grow randomly. In 
the region 1 < ^ < 2.1, the attenuation vector increases with the increase in z. 
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face Plot of Attenuation Vecto 




Surface Plot of Refractive index 



Surface Plot of 





Figure 6: The waves move towards the event horizon. Most of the region 
admits the properties of anomalous dispersion 



The variation of refractive index with respect to uj is negative in most of the 
region which shows anomalous dispersion of waves. There are some points in 
the region where ^ > and dispersion is normal. Figure 7 indicates that the 
attenuation vector takes negative values. The propagation vector decreases 
with the increase in angular frequency and z. The attenuation vector takes 
random values with the increase in angular frequency. The refractive index 
is greater than one and its variation with respect to angular frequency admit 
positive values at some points which indicates normal dispersion of waves. 
Figure 8 shows that the attenuation vector is positive throughout the region. 
The propagation vector decreases with the increase in angular frequency and 
z. The attenuation vector randomly increases and decreases with the increase 
in angular frequency. It shows that the waves randomly damp and grow 
with the increase in angular frequency. The refractive index is greater than 
one and ^ > in most of the region which results that the waves disperse 
normally. Figure 9 shows that the attenuation vector is positive for the region 
1 < 2; < 5, 10~^ < C(j < 10. In a small region near the event horizon, the 
attenuation vector decreases with the increase in z which indicates that the 
waves damp as they move towards the event horizon. The refractive index 
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Figure 7: The wave are directed towards the event horizon. The dispersion 
is found to be normal and anomalous randomly 

Surface Plot of Propagation Vector 




Figure 8: The waves move towards the event horizon. Most of the region 
admits normal dispersion 
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Surface Plot of Propagation Vector 
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Surface Plot of Attenuation Vector 
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Figure 9: The waves are directed towards the event horizon. Some of the 
region admits normal dispersion 

is greater than one and its variation with respect to angular frequency is 
greater than zero at random points which admit normal dispersion of waves. 

In Figures 6-9, the propagation vector is negative in the whole region. 
Moreover, Figures 6 and 9 admit a small region near the event horizon where 
the refractive index increases with the decrease in z. 

7 Non- Rotating Background with Isothermal 
Plasma 

For the isothermal plasma model, the Fourier analyzed perturbed GRMHD 
equations (3.3.1)-(3.3.5) of ^26j are given as follows 






(7.1) 
(7.2) 



a 

Lkc5 = 0, I 
ci{—pLuj + ikpau — (uap)' — au^'j'^pu'} + C2{—pLuJ + ikpau 
+ {uap)' + a;M^7^pu'} + c^^p + p){a{l + '-y'^u'^)tk — a(l — 27^u^) 
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u' 



xil + Yu^) i0JYu} = Q, (7.3) 

u 

cip-f^{a, + uu\l + 7V)} + C2{p-f^{a, + uu\l + 7V)} 



+Lkp + p'} + cs{p + p)7^[(l + l^u^){-— + iku) + {u'{l + 7^2' 

a 

x(l + 47V) + 272MaJ] = 0, (7.4) 

ciP7 I \- u[az + 7 )| + C2p{ (7 - 1) + 7 u[a^ 

a a 

+7^-u-u')} + c-iip + p)7^{ 7^u + j'^u'^ik + + a^) 

a 

x(l + 27V)} = 0. (7.5) 

Equations (17. ip and (17. 2p show that C5 = which means that there are no 
perturbations in magnetic field. 

7.1 Numerical Solutions 

In order to find the numerical solutions we assume the same time lapse and 
z-component of fiuid velocity as in Section 4.1. When we substitute these 
values with the assumption of stiff fiuid, i.e., p = p, the mass conservation 
law in three dimensions gives P = = P- The complex dispersion relation 
is of the form 

Ai{z)k'^ + A2{z, uj)k + A^{z, uj) + L{Ai{z)k^ + A^{z, uj)k'^ 
+A^{z,uj)k + Aj{z,uj)} = (7.6) 

This gives three complex values of k shown in Figures 10-12. 

Figure 10 indicates that the propagation decreases while the attenuation 
vector increases with the increase in angular frequency. The refractive index 
is greater than one and its variation with respect to uj is greater than zero 
in the region 0.5 < z < 2. Rest of the region admits points of normal as 
well as anomalous dispersion. Figure 11 shows that the propagation vector 
decreases with the increase in angular frequency. The attenuation vector is 
also negative and it decreases with the increase in angular frequency and z. 
It shows that the waves grow with the increase in angular frequency and z. 
The refractive index is greater than one and its variation with respect to 
UJ is negative which indicates that the region admits anomalous dispersion. 
Figure 12 shows that the attenuation vector admits negative values. The 
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Surface Plot of Propagation Vector Surface Plot of Attenuation Vector 




Figure 10: The waves move towards the event horizon. A small region near 
the event horizon admits normal dispersion 



Surface Plot of Propagation Vector Surface Plot of Attenuation Vector 




Figure 11: The waves are directed towards the event horizon. Dispersion is 
found to be anomalous 
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Surface Plot of Propagation Vector 




Surface Plot of Attenuation Vector 
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Figure 12: The waves move toward the event horizon. The dispersion is 
normal in a small region near the event horizon 

propagation vector decreases while the attenuation vector shows random in- 
crease and decrease with the increase in angular frequency. The attenuation 
vector decreases near the event horizon which indicates that the waves grow 
near the event horizon. The refractive index is greater than one. The change 
in refractive index with respect to angular frequency is greater than zero in 
the region 0.5 < z < 3, 0<u;<10 which indicates normal dispersion. 

The propagation vector takes negative values for Figures 10-12. In a 
small region near the event horizon, the refractive index increases with the 
increase in z for Figures 11-12. 



8 Rotating Non-magnetized Background with 
Isothermal Plasma 

The Fourier analyzed GRHD equations for isothermal plasma are given by 
Eqs.(3.4.1)-(3.4.4) of [26]. 





Ci{(— to; + iaku)p — a'up — au'p — aup' — a'j^up{VV' + uu')} 
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+C2{{—LU! + Laku)p + a'up + au'p + aup' + a^'^upiVV' + uu')} 

u' 

+C3(/0 + p)[-6Cc;7\ + tfcQ;(l + 7^^) - a{(l - 27\^)(1 +7V) — 

-2-i\i'^VV'}] + C4(p + p)[^'^V{-LUJ + iA;^^) + a^'^u{{l + 27V2)V' 
+272?7,\/?/}] = 0, (8.1) 
cip7\{(l + 7V2)V^' + ^"^uVu'} + C2P7^m{(1 + -i'^V^)V' + 7^,1/?/} 



+C3(p + p) 



f ( ^ + ) + {(1 + 27V2)(1 + 27^x^2) _ ^•^v^}v' 



a 



+2-i'^uVu'{l + 272'u2)] + Ci{p + p)7^ 



q; 



+7^ii{(l + 47^1*^)^ + 4(1 + 7V^)yy'}] = 0, (8.2) 
CiP7^{az + W(l + 7^w^) + 7^w^yy} + C2[p7^{a^ + W(l + 7^^^) 



+YuVV'} + p' + i/cp] + C3(p + p)7 



(I + 7V) 



iku 



+u'{l + 7\^)(1 + A-fW) + 2V{(1 + 27^ii2)yy' + a J] 



+C4(p + p)7^ 



— io; 



a 



+ iku]+ u'V'{l + 47^1x2^ 



+2V^{(1 + 272^2)^^' + a,}] = 0, 



(8.3) 



cip7 



—to; 
a 



+ 4a^ + 7'(yy' + W)} 



+ C2P 



'''^(7^-1)+ 



a 



+7^(yy + W)}] + 03(^ + ^)7^ 17^'^ (^-^ + ''^w^ + {az + 7^'"^i') 

(1 + + 7Vl/'(l + A'^W)] + C4(p + p)7^ I7V (^-^ + 

+7^1^^' (1 + 47V2) + 272-^^(0^ + 27^^) } = 0. (8.4) 



8.1 Numerical Solutions 

Using the same assumptions as given in Section 5.1, we obtain a complex 
dispersion relation of the form 

Ai{z)k'^ + A2{z, uo)k^ + Ay,{z, uj)k'^ + Ai{z, uo)k + ^5(2, uo) 
+i{AQ{z)k^ + A-ji^z, u)k'^ + A^{z, uj)k + Ag{z, cj)} = 0. (8.5) 

The four complex values of k are represented in Figures 13-16. 
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Surface Plot of Propagation Vector Surface Plot of Attenuation Vector 




Figure 13: The waves move towards the event horizon. A small region near 
the event horizon admits anomalous dispersion 

Figure 13 shows that the propagation vector takes negative values. The 
propagation vector decreases with the increase in angular frequency and z. 
The attenuation vector increases with the increase in angular frequency. The 
refractive index is greater than one and its variation with respect to lo is less 
than zero in the region 0.75 < z < 4.5, < < 10 which leads to anomalous 
dispersion. Random points of normal and anomalous dispersion are found 
otherwise. Figure 14 indicates that the propagation and attenuation vectors 
decrease with the increase in z and lo. It shows that the waves grow with 
the increase in angular frequency and z. The refractive index is greater than 
one and its variation with respect to lo is positive in the region 0.75 < z < 
4, < < 10 which indicates normal dispersion of waves. Random points of 
normal and anomalous dispersion lies in rest of the region. Figure 15 shows 
that the propagation vector decreases with the increase in angular frequency 
and z. The attenuation vector decreases as z decreases which indicates that 
the waves grow as they move towards the event horizon. The refractive index 
is greater than one. The quantity ^ > at random points which leads to 
normal dispersion of waves at those points. In Figure 16, the attenuation 
vector decreases as z decreases in the region except for a; = which shows the 
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Surface Plot of Propagation Vector 



Surface Plot of Attenuation Vector 





Figure 14: The waves are directed towards the event horizon. Normal dis- 
persion of waves is found near the event horizon 




Surface Plot of Refracrive Index Surface Plot of — — 




Figure 15: The waves are directed towards the event horizon. Dispersion is 
found to be normal at random points 
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Figure 16: The waves move towards the event horizon. Dispersion is found 
to be normal at random points 

growth of waves as they move towards the event horizon. The refractive index 
is greater than one and its variation with respect to u is greater than zero 
at random points which shows that the waves disperse normally at random 
points. 

In Figures 13-16, the propagation vector admits negative values which 
indicates that the waves move towards the event horizon. The refractive 
index increases in a small region near the event horizon for all these figures. 

9 Rotating Magnetized Background with 
Isothermal Plasma 

This is the most general and physical case considered to investigate the wave 
properties, i.e., we take rotating isothermal plasma which is obviously highly 
magnetized. The respective Fourier analyzed perturbed GRMHD equations 
are given by Eqs.(3.5.1)-(3.5.7) of 

Ci{a' + ika) — CaKaA)' + ikaX] + c^{aV)' — CQ{{au)' — luj 
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+Lkua} = 0, (9.1) 

C5(^^ + ika) ^0, (9.2) 

a 

C5ik = 0, (9.3) 
Ci{(— /.cu + iaku)p — a'up — au'p — aup — a'y'^up{VV' + uu')} 
+C2{(— to; + iaku)p + a'up + au'p + aup' + a'y'^up{VV' + uu')} 

+C2,{p + p)\-iu^'^u + iA;Q;(l + 7V) - q;{(1 - 27V) (1 + 7^^^) — 
-27^iiVy'}] + C4(p + p)[7V(-ta; + ikau) 

+a^'^u{{l + 27V2)l^' + 2^'^uVu'}] = 0, (9.4) 
Ci7V^i{(l + -f^V^)V' + 'y^Vu'} + C27V{(1 + -f^V^)V' + ^\Vu'} 

XB"^ luj XB^ 
+C3[-{(p + p)l^uV - -—}— + {{p + p)l^uV + -—}iku 

47r a 47r 

+ (P + P)7'{(1 + 27'7i2)(l + 27V2) - 7 + 2(p + p)-f\l 

+2j^u^)uVu' + P^iaX)'] + C4[-{(p + p)j'{l + 7V2) + 

47ra 47r a 

+{(P + P)7'(l + 7^2) - + (p + p)7'«{(l + 472^2)^^' 

+AVV'il + 7V^)} - - ^C6{(1 - u')Lk + -(1 - u') 

Ana An a 

-uu'} = 0, (9.5) 
cip-f^{a, + (1 + 7^m^)mm + 7^MVy'} + C2[p7^{a^ + (1 + -f^u^)uu' 

+-f\W} + ikp + p'] + C3[-{(p + p)7'(l + 7V) + 

An a 

+{{P + P)7'(l + 7V) - + {(P + + 7'«') 

(1 + A^^u^) + 2u^\a, + (1 + 2^^u')VV')} - ^(aA)'}] 

+C4[-((p + p)7^'«^- ^)— + {(P + P)7^'«^+ ^I'-^w 

An a An 

+{(P + p)7'{«V'(l + 47V2) + 2\^(a, + (1 + 2^\^)uu')} + h^^}] 

Ana 

+— C6[A(1 - u'')Lk + A-(l - u^) - Xuu' + = 0, (9.6) 

An a a 

CiP7^[-— + u{a^ + -t'^{VV' + W)}] + C3[- — {(p + p)7^M 
a a 
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+ -^{uX -V)} + Lk{{p + p)7V + -—} + (p + ph^{{a, 

+27W)(1 + 2^\') + 7VV'(1 + A^\')} + ^(3Aa, + A')] 

+C4[- — {(p + ph'V -^(u\-V)} + ik{{p + p)^'Vu - ^} 
a An An 

+ {p + ph\{V'{l + + 2V{a, + 27^^)} - -—aA] 

An 

+C6 — [ {A + 2m(uA - V)] + LkuX + 3aJ2Au - + 2m'A 

An a 

+mA' - V] + C2p[- — (7' - 1) + M7'{«. + 1^{VV' + W)}] = 0. (9.7) 
a 

Equations (19.21) and (19. 3p give C5 = which shows that the 2;-component of 
magnetic field is not affected by gravity. 

9.1 Numerical Solutions 

We assume the same values of time lapse, x and ^-components of the fiuid 
velocity, A and B as given in section 6.1. Using these values with the assump- 
tion p = p, the mass conservation law in three dimensions gives P = = P- 
We obtain a complex dispersion relation quintic in k, i.e., of the form 

Ai{z)k^ + ^2(2:, uj)k^ + A^{z, uj)k^ + ^4(2:, uj)k + A^{z, uj) 
+L{AQ{z)k^ + Ai{z, uj)k^ + A^{z, uj)k^ + Ag{z, uj)k'^ + ^10(2, uj)k 
+An(2,^)} = 0, (9.1) 

which can not be solved to get the exact solutions. We have found the 
numerical solutions with the help of software Mathematica. We obtain five 
complex values of k shown in Figures 17-21. 

Figure 17 shows that the propagation vector decreases with the increase in 
angular frequency. In a small region near the event horizon, the attenuation 
vector increases as z decreases. The refractive index is greater than one 
and its variation with respect to u is positive at random points which shows 
normal dispersion of the waves. Figure 18 shows that the propagation vector 
decreases with the increase in angular frequency. The attenuation vector 
randomly increases and decreases in the region l<z<2.2, 0<ci;<10. 
The refractive index is greater than one and its variation with respect to uj is 
positive at random points which shows that the waves disperse normally at 
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Figure 17: The waves are directed towards the event horizon. Normal dis- 
persion is found at random points 





Figure 18: The waves move towards the event horizon. The waves disperse 
normally at random points 
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Surface Plot of Propagation Vector Surface Plot of Attenuation Vector 




Figure 19: The waves move towards the event horizon. Dispersion is found 
to be normal at random points 



Surface Plot of Propagation Vector Surface Plot of Attenuation Vector 




Surface Plot of Refractive Index Surface Plot of 

did 




Figure 20: The waves move towards the event horizon. Normal dispersion is 
found in a small region near the event horizon 
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Surface Plot of Propagation Vector 





Figure 21: The waves are directed towards the event horizon. Normal dis- 
persion is found at random points 

those points. Figure 19 indicates that both the propagation and attenuation 
vectors are negative and decrease with the increase in angular frequency. 
The refractive index is greater than one and its variation with respect to lo 
is positive at random points which shows that the waves disperse normally 
at those points. Figure 20 indicates that both the propagation and the 
attenuation vectors decrease with the increase in angular frequency. The 
refractive index is greater than one and its variation with respect to uj is 
greater than zero in the region 1 < z < 1.2, 0<ti;<10 which shows normal 
dispersion. Random points of normal and anomalous dispersion are found 
otherwise. Figure 21 shows that the propagation vector decreases with the 
increase in angular frequency. The attenuation vector increases and decreases 
randomly in the region l<z<2.2, 0<a;<10. The refractive index is 
greater than one and its variation with respect to uj is positive near the event 
horizon which indicates normal dispersion of waves. 

In Figures 17-21, the propagation vector admits negative values indicat- 
ing the direction of waves towards the event horizon. Furthermore, in all 
these figures, a small region near the event horizon admits the increase in 
refractive index with the decrease in z. 
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10 Outlook 



We have discussed the cold and isothermal plasma wave properties of the 
Schwarzschild black hole magnetosphere in 3+1 formalism. 

For this purpose, we have used the planar analogue (given by Eq. fl2.1l) ) 
due to the following reasons: 

• It is difficult to study the MHD stream equation and the behavior of 
a perturbed magnetosphere in the Schwarzschild spacetime. We have 
chosen a different way and altered the spacetime (i.e., we have taken 
the planar analogue) which simplifies the analysis of magnetosphere. 
It has been done while preserving the following key features of the 
Schwarzschild metric: 

1. The behavior of the Schwarzschild time coordinate as r — ^ 2M is 
analogous to the behavior of Schwarzschild planar analogue time 
coordinate as z ^ 0. 

2. The planar spacetime and its stationary MHD magnetosphere 
have nice computational features that all aspects of the magneto- 
sphere should become asymptotically z independent at large z. 

• The planar analogue is not empty which is not important for our anal- 
ysis because it serves merely as a testbed for studying various aspects 
of the interactions of gravity with plasma [20]. Therefore, we are free 
to assume that there is no direct non-gravitational interaction between 
the plasma and materials whose stress-energy produce the spacetime 
curvature. 

The plasma is assumed to be living in this planar analogue for which the 
dispersion relations are obtained from the Fourier analyzed perturbed 3+1 
GRMHD equations for the non-rotating, rotating non-magnetized and rotat- 
ing magnetized backgrounds. These relations give the value of z-component 
of the wave vector which yields the propagation and attenuation vectors, 
the refractive index and its change with respect to angular frequency. The 
graphs are obtained to have a clear viewpoint. A summary of the results is 
given below. 

For the cold plasma living in the vicinity of Schwarzschild regime, nor- 
mal dispersion of waves is found in the Figures 1, 3 and 8 which indicate 
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that the waves can pass through the medium. Figures 2, 4, 6 and 9 rep- 
resent that the waves disperse anomalously in most of the region. Random 
dispersion of waves is found in the Figures 5 and 7. Figure 3 shows that the 
waves damp in a small region near the event horizon. In Figures 4 and 5, 
the waves grow near the event horizon in a small region. In Figure 9, the 
attenuation vector increases with the decrease in 2; in a small region near the 
event horizon. In Figures 1-9, the waves damp and grow randomly with the 
increase in angular frequency. 

For the isothermal plasma surrounding the Schwarzschild black hole, 
Figure 14 admits normal dispersion near the event horizon. The waves dis- 
perse anomalously in Figure 11. In Figure 10, the attenuation vector in- 
creases whereas it decreases for Figure 11 with the increase in angular fre- 
quency. It shows that the waves damp in Figure 10 and grow for Figure 11 
with the increase in angular frequency. Figure 12 shows that the wave grow 
and damp randomly with the increase in angular frequency. In Figures 15 
and 16, the waves grow in a small region near the event horizon. In Figure 
17, the waves damp whereas they grow in Figures 18 and 21 with the de- 
crease in z in small regions near the event horizon. In Figures 19 and 20, 
the attenuation vector decreases with the increase in angular frequency. 

For the cold plasma living in the non-rotating Schwarzschild background, 
there exists a case where the waves disperse normally (Figure 1) whereas for 
the isothermal plasma, normal dispersion lies at random points. In the rotat- 
ing non-magnetized background, the cold plasma shows normal dispersion in 
Figure 3 whereas the isothermal plasma admits normal dispersion in a small 
region near the event horizon in Figure 14. In the rotating magnetized back- 
ground. Figure 8 indicates that the cold plasma admits normal dispersion in 
most of the region whereas the isothermal plasma admits normal dispersion 
at random points. Thus we can conclude that the pressure ceases the normal 
dispersion of waves. 

In most of the cases, the refractive index increases as the waves move to- 
wards the event horizon in a small region near the event horizon. This shows 
that the refraction of waves increases as they move towards the event hori- 
zon. For all the backgrounds of cold and isothermal plasma, the propagation 
vector takes negative values which shows that the waves move towards the 
event horizon. 
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